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Einstein Condensate qualitatively illustrate the weak and strong regimes respectively. 
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1. Introduction 

Decoherence effects are now believed to be the essential ingredient which destroys most of 
the counterintuitive aspects of quantum mechanics. Such effects are at the same time an 
academic tool to the understanding of the classical limit of quantum mechanics as well as 
an important ingredient in the area of quantum computation. The dynamics of quantum 
open systems has been therefore extensively studied [1]. Of particular importance in this 
context is the Born-Markov approximation which leads to master equations of various 
kinds |2], whose validity is limited by the weak coupling approximation. The strong 
coupling regime however has been less explored. 

It is the purpose of the present contribution to shed some light onto the weak- 
strong coupling transition in the context of two linearly interacting systems under the 
influence of a phase fluctuating coupling. The model in spite of its schematic character 
has been shown in several instances and different areas to reflect and adequately 
describe experimental results. Examples are the description of exciton-polariton damped 
oscillations [3l H]; predictions for the behavior of oscillations of two coupled modes in 
the context of microwave cavities [H [6]; the self-trapping phenomenon in the tunneling 
process of a Bose-Einstein condensate (BEG) [7l|8l|9llini[IIl[12l[13l[Tl]. 

We will show that in the weak coupling regime the usual master equation results are 
recovered and the usual phenomenological damping constant is derived as a function of 
the model parameters. The strong coupling limit however leads to a completely different 
scenario, the "freezing" of the dynamics. We will illustrate these effects in the context of 
exciton-polariton oscillations and the self-trapping of a BEG in a devised laser potential. 

2. The Model 

In this section we give a detailed derivation of the stochastic time evolution of the 
following system 



where a"'' (a) and b^{b) are creation (annihilation) operators. Ua and Ub are the frequencies 
of the a-mode and b-mode respectively. Since that the operators a and b were considered 
as bosons they will obey the commutation relation for bosons [a, a^] = 1 and [b, b^ = 1. 
In the third term, g{t) stands for the interaction strength between the modes and here 
is assumed to be time dependent in the sence that \git)\ is constant but its phase is a 
stochastic variable. The Hamiltonian ([T]) can be written in matrix form 



H = huJaO^a + huJbb'^b + h[g{t)a)b + g*{t)b^a] 



(1) 




(2) 



The time evolution operator for the system is given by 




cos(l^lt) ^sin(|(?|t) 
^sin(|(?|t) cos{\g\t) 



) 



(3) 
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where Hint is the interaction part of the Hamiltonian H. 

We will specify the noise by defining a stochastic process for (j){t). In our model we 
assume that 

g{t) = goexp[i(j){t)], (4) 

where go is the non-stochastic amplitude while the phase is treated as a stochastic 
variable. Here, we will consider random phase telegraph noise where itself fluctuates 
in the manner of jumps. In particular, the phase fluctuations were describe by a 
Wierner-Levy (phase diffusion) process and the amplitude fluctuations by a colored 
gaussian noise. An alternative model which represents noise by means of discrete jump 
processes was first introduced into quantum optics by Burshtein and Oseledchik jl5| . 
A simple example of such a jump process is the two-state random telegraph. These 
models are very convenient and elegant to study the noise of the electromagnetic atom- 
field interaction in a non-perturbative manner. The random telegraph models, whether 
associated with phase, frequency or amplitude fluctuations lead to an equation for 
average responses in exact algebraic form. The model of random telegraph (jump-type) 
noise is physically very sound to describe the noise arising from electromagnetic field 
fluctuation or from collisions of various kinds or from other external sources. Indeed, 
that model including the effects of stochastic phase and/or in amplitude has been 
explicitly solved for the case of the James-Cumming model (JCM) by A. Joshi |16j . 
The fluctuations are modeled by the random telegraph process and an equation for 
the density operator averaged over the fluctuations is obtained. The solution of these 
equations was used to study the decoherence effects in the dynamics of the system. A. 
Joshi's work was treated in a pedagogical form by E. A. Ospina [I7j. We assume further 
that the change in (j)(t) occurs instantaneously jump wise and the jumps are separated 
by mean time intervals of the order tq in which 0(t) = constant, as shown in Fig. ([T]). 
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Figure 1. Schematic representation of the phase distribution of coupling between the 
modes. 

There are two stochastic variables: the time interval r between one and the next 
jump and the value of the phase constant (p in each of these intervals. The variable 
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T = ti — ti^i follow the probability distribution 

dQ{T) = -e-"/"«c/r. (5) 

^0 

with to = 0. The above distribution specifies the probability of duration of each such 
jump interval and has intervals mean duration tq. We consider only the case in which 
the phases are uncorrelated. The probability distribution to phase is given by 

dq{<P) = (6) 

with mean value (0) = vr. So, at any instant, the probability of finding a given 
remains the same and equation to dQ{t), and there is no limitation on the form of this 
distribution. In other words, (j){t) is undergoing random continuous change of Markov 
type. 

The dynamics of the system is given by the unitary transformation U (0, t, t ) such 

that 

p{t-(P) = u{4>,t,t)p{t)u-\4>,t,t). (7) 

At the end of each (zth) interval we find the density matrix p(t) which is the initial 
condition for the next matrix, so, if in the interval (0,t) there are k jumps in 0, then 

p{t; ti, ■ ■ ■ , tfc , 00, ■ ■ ■ , 0fc) = (0fc; t, tk)U{(f)k^i; tk, tk-i) ■ ■ ■ 

X f/(0i; t2, ti)f/(0o; ti, O)p(O)f/-i(0o; ti, O)f/-^(0i; ts, ti) ■ ■ ■ 

X U-\(l)k-utk,tk-i)U-\(j)k]t,tk). (8) 

The above expression is of multiplicative nature and hence it is quite easy to average 
over. The probability (in the interval (0,t)) that k changes of have actually occurred 
at successive instants ti,t2, ■ ■ ■ ,tk and that a certain sequence of 0i, 02, ■ ■ ■ , 0^ (where 
(pi = (p(ti)) was realized between them is obviously equal to 

dP{tu t2, ■ ■ ■ , tk, 01, 02, ■ ■ ■ , t) = ^e~*/^" ( J] rfti) ( J] dg(0,)) . (9) 

'^0 i=l i=l 

The average density operator can thus be written as 



k=0 



P{t) = ^ I I dP{ti, t2, ■ ■ ■ , 4; 01, 02, ■ ■ ■ , 0fc, t)p{t; ti, ■ ■ ■ , tfc, 00, ■ ■ ■ , 0fc). 

(10) 



Rewrite ( !T0|) with use of ([9]) we get 



-p{t)e 

k=0 



y2~k dtk-^i--- I dti / dq{(f)k) 

j,_n '^0 ^0 ^0 Jo J 



X y rfg(0o)p(t;ti,---,tfc,0o,---,0fc). (11) 

Note that the term with k = (when does not change at all in the interval (0, if:)) will 
not contain integrals with respect to time and thus is given by 

dg(0o)p(t,0o). (12) 
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Now using the recurrence relation ([7]), we can multiply both sides of equation fill I) 
from the left (right) by U~^{(j)] T,t) (U^cj)] T,t)) respectively and also by dq{(f))dt/TQ, then 
integrate with respect to time from to r and eliminate the entire series using equation 
( ITTj) . After some simplifications it is easy to show that 

p(r)e^/^»= [ dq{<Po)U{<l)o;r,0)p{0)U~\<Po;T,0) 



+ - [ dte'/^'' [ dq{(t))U{(P]T,t)p{t)U-\(P;T,t). (13) 

To Jo J 

Now we will rewrite the equation (|T3|l above in term of matrix elements of operators 
p(r) and U (0; r, t) 

k,l 

+ 1 r dte"^' [ dq{<l))y^Ui<l)-T,t),kP{t)MU-\<l);T,t)i^.{U) 

Define the conjunct of matrix {G*™(r, t)} with elements give by 

G'"(r,t);,= j dq{(P)UMr,t)Urr^{(P;T,t), (15) 

so 

/ k 

+ 1 Tdte^Z-J^ [5^G^'"(r,t),fep(t),/ 
-'o Ik 

= ^[G^-(r,0)p(0)], 
I 

+ 1 TdteVn, \2prn^r,t)p{t)]u. (16) 

roJo V 
Using the trace definition TrA = S^^ 

p(r)i^ = e--/-«Tr[G^-(r,0)p(0)] 

+ -/ rfte-("-*)/"°Tr[G*'"(r,t)p(t)l. (17) 
^0 Jo J 

This is the statistical average over the random variable 4>{t). In order to determine 

the dynamical evolution of the system one has to determine G. The problem is now 

simplified because we have to deal with interval in which (p (or alternatively g) is constant 

and the change in p is perfectly regular. Thus, knowing G, we can find the average 

variation of the system during the relaxation process [15]. To determine G we will use 

the time evolution operator U{t) given by ([3]) with g(t) defined in (jlj) 

r/m = I cos(|^|t) e^<^sin(|^|t) . 

I _e-^4>sm{\g\t) cos {\g\t) ' ' ^ ' 



Stochastic coupling in two modes systems ... 6 

The elements to be averaged in the calculations of G are those containing the factors 
e^*''^. Since the phases are equally probable most of the terms vanish after averaging. 
The remaining (relevant for our purposes) non-vanishing elements of 

G\\ = G'^ = cos''[{g,{r-t)] 
sin^ [{g^{T-t)] 



are 



rill _ r^22 



^22 
ril2 



11 
-y21 



(19) 



G?^ = cos2[(^o(r-t)]. 



We have thus 



G'\r,t) 



G'\r,t) 



G''{T,t) 



cos^ [(c/o(r - 1)] 











sin^ [(^o(r - t)] 





cos2 [((7o(r - t)] 
sm^igoir-t)] 



[G'\r,t)y 



(20) 







cos2 [{goir - t)] 

and using equations ( IT71) and (1201) we obtain 

p(r)ne^/^° = p(0)ii + [p(0)22 - p(0)ii] sin^ (gor) 

+ 1 /" dte*/"°{p(t)ii + [p(t)22 - p(t)ii] sin' h(r 



(21) 



p(r)22e"/"° = p(0)22 + [p(0)ii - p(0)22] sin' (g 



'or 



+ - rfte*/"°{p(t)22 + [p(t)ii - p(t)22] sin' [goi 
^0 Jo 

3. Dynamics of the average number of a mode 

The average bosons number is given by 

na{t) = Tr p{t)a)a 
and can be evaluated using equations fl^T]) : 

Mr))e^'^' = {nam [l - 2 sin' ((^qt)] + (N) sin' (gor) 

3*/""(M^)|l - 2sin' [go{T - t)]]dt 



(22) 



+ 



^0 Jo 

(N) 



(23) 



e*/^« sin' [^o(^ - t)]dt, 



To Jo 

where (A^) = {na{0)) + {nb{0)) is the initial excitation number with {na{0)) and {nf,{0)) 
being the average excitations in each mode. This equation describes the relaxation of 
the intensity of the mode a and can be solved using the Laplace transform. To solve 
(!23|) we will define 



fit) = (n„(t))e*/-°, 
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g{t) = 1-2 sin^ (^fot) = cos {2gQt), 
hit) = e*/"°, 
j{t) = sin^ (got). 

inserting the functions defined above into (123!) we obtain 



fir) 



M0))9ir) + ^ , 
^0 Jo 



+ {N)j{r) + 



(AO 

To 



f{t)g{T - t)dt 
h{t)i{T - t)dt. 



(24) 



Now applying the Laplace transform on both sides of (|24l) we obtain 

f{s) = {riamgis) + -f{s)g{s) 

To 



+ {N)j{s) + ^-^h{s)j{s) 

To 



(25) 



{nM)9is) + {N)jis) + mis)jis) 



1 - Tjis) 



Calculating the Laplace transform of the functions g{s), h{s) and j{s) and substituting 
it into (l25l) we obtain 

^2 



fis) 



where we defined Q 



K(o))s 



{N)i2go 



2ro J 



+ 

+fi2 2[s- 



1 



2ro J 



V 



+ ^2 



(26) 



i^go 



(2roy 



■ . The inverse Laplace transform of Eq. fl26p yields 



the following expression the time evolution of the relaxation as 

sin (Qt) ' 



Mt)) 



+ 



K(0))e-*/2"° cos (fit) + 



2rofi 
*/2-o _ cos (fit) 



sm 



2rofi /■ 



(27) 



In the limit when the average time tq between phase jumps is large as compared with the 
oscillation period one obtains for goTo — )■ oo a pure oscillatory regime with frequency 
2gQ. As the average time between frequency jumps increases one obtains an envelop 
limiting the oscillation amplitudes and the oscillation frequency is only slightly altered. 
However in the limit where the average time between jumps decreases as compared 
with the coupling g^ the oscillation ceases at goTQ = 0.25 and an over damped limit 
sets in. Figures ([2]) and (EI) illustrate the weak coupling regime {goTo ^ 0.25) and 
strong coupling regime {goTo > 0.25) respectively. Figure @ illustrates the "freezing" 
of the dynamics as goTo decreases beyond the limit goT^ = 0.25. The fluctuating interact 
records information of the system. The information transfer plays the role of unobserved 
detection process [H]. In our case small values of g^To would imply according to this 
reasoning that the system is being "measured" with increasing frequency, freezing out 
as a Zeno like effect. 
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(^^(t)) = J e*/2-o _ cos m] - (2^ 



4. Applications 

4.I. The weak coupling regime: Exciton-polariton oscillations 

Let us consider the weak coupling regime (WCR), where (70 To > 0.25, and consequently 
is a real number. Note that, as tq — )■ 00, the expression for {na{t)) reduces to the 
usual result without fluctuations. The effects of phase fluctuation in the intensity of the 
mode a , e.g, for an initial number state |v9(0))f, = \Ni,) and the mode a in vacuum state 
are given, according to (ITTI) . by 

The result above shows that the intensity of the mode a contains two parts: (1) Rabi 
oscillations with frequency Vl] (2) a comparatively slow- varying part e~*/^^°. To see this 
more clearly, we plot equation fl28|) in Figs. ([2]) and ([3]) for A*";, = 2. It is clear that 
the damping of Rabi oscillations is more pronounced when the mean time interval tq 
between phase jumps become shorter and shorter until g^TQ = 0.25 (see Figures ([2]) 
and ([3])). In other words the decoherence mechanism is faster for shorter-phase jump 
intervals. 

In this section we compare the results obtained in Ref. [4] which describe the 
exciton-polariton oscillations in the weak coupling regime with a phenomenological 
coupling constant. In a real cavity, the modes of exciton and photon are coupled to 
a continuum of modes which leads to dissipation. The coupling can be scattering of 
phonons in the case of exciton or cavity damping in the case of photon. In both cases, 
the result is to dampen the mode of interest. The result obtained in Ref. |1] can be read 
of Eq. (l28|l where the damping coefficient corresponds to Tq = = (7p + 7ex)~^, where 
7e2: and 7p are exciton and photon damping from the reservoirs. Using the values of 7 
adopted in Ref. |1] we may estimate the order of magnitude of tq as 10~^°s. The model 
has been shown to reproduce experimental results [HI [20] . In this situation equation 
( l28i) can be written as 



M)) = Y [1 - e-^*/' cos {2got)\ . (29) 

The results obtained here describe the decoherence process in the system. However, 
the fluctuations introduce a flnite width in the transmission spectrum even in a lossless 
cavity. On the other hand, in recent work, S. Schneider et al., [2T] also included fluctu- 
ation in intensity and phase in the exciting laser pulse to explain effects of decoherence 
for single trapped ion. In Schneider's model the intensity and phase fluctuations define 
a stochastic Schrodinger equation in the Ito formalism [22], or more appropriately a 
stochastic Liouville-von Neumann equation. The results are in good qualitative agree- 
ment with recent ion experiments 
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2,00 




Figure 2. Light intensity as a function of dimensionless time T = got for the case 
that the excitons are initially in a number state N ~ 2 for qoTq — ?> oo (dotted line), 
goTQ = 100 (dashed line) and g^TQ = 10 (solid line). The intensity is in arbitrary units. 




Figure 3. Light intensity as a function of dimensionless time T = got for the case 
that the excitons are initially in a number state Nf, ~ 2 for g^T^ = 1.0 (dotted line), 
goTo — 0.50 (dashed line) and 5oto = 0.25 (solid line). The intensity is in arbitrary 
units. 



4-.2. An alternative self-trapping mechanism 



Now we analyze the strong coupling regime (SCR), where qqTq < 0.25. In this case Q is 
purely imaginary and Eq. fl27|) can be written as 



Mt)) 



JV^e-i/2ro 



,t/2T0 



cosh{\n\t) 



smh{\n\t) 
2to\Q\ 



(30) 
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The SCR may be investigated by looking at the intensity of the mode a. As observed 
above when qqTq becomes shorter and shorter as compared to 0.25, the fluctuation 
effects are larger, and the fluctuations prevail over the oscillation between mode a and 
mode b. In this case the SCR modifles the picture. The inhibition of the transition of 
excitations between the modes is induced by the fluctuations in the coupling. This 
can be interpreted as an environment induced "quantum Zeno-like effect (QZLE)" 
O [211 [25], [26], [271 [28]. In the regime QoTq <^ 0.25 the interaction between mode a and 
mode b is not able to absorb or release energy and therefore stay put. The fluctuations 
in the interaction strength between the mode a and mode b inhibits the excitation of 
mode a (in Fig. (jlj), we exemplify this effect). When g^ro — )■ 0, {riait)) — )■ 0, when 
(^a(O)) = 0; the dynamics is frozen. In Ref.[2| a self-trapping mechanism of BEC in a 
laser potential has been reported. Two explanations have been given. Firstly the one 
using a nonlinear Gross-Pitayesty equation P, |TT] and the other a schematic many body 
system [H]. In the present contribution one might view modes a and b as the two sides 
of the well and the self-trapping mechanism as the freezing out of the dynamics due to 
uncontrollable fluctuations in the experiment. 
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Figure 4. Light intensity as a function of dimensionless time T — g^t for tlie case 
that the excitons are initially in a number state — 2 for qoTq = 0.25 (solid line), 
goTo = 0.01 (dashed line), goTo = 0.001 (dashed dotted line) and 50^0 = 0.0001 (dotted 
line). The intensity is in arbitrary units. 



5. Conclusion 

We studied a system of two linearly coupled oscillators and the effects of a phase 
fluctuating coupling. The model can be solved analytically and displays the weak- 
strong coupling transition. We show that this transition is a function of a dimensionless 
parameter goTo and occurs at goTo = 0.25. In the weak coupling regime we provide for 
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an analytical expression for the damping parameter and compare with that of Ref . |1] , 
in the context of exciton-polariton oscillations. The strong coupling regime leads to a 
"freezing" of the dynamics and may qualitatively provide for yet a third explanation for 
the self-trapping phenomenon in BEC (the first two given in Refs. [9l[TT] and [l4j). 
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